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ABSTRACT 


This  study  is  concerned  witn  the  stabilization  in  attitude  of  an  earth 
satellite.  In  particular  two  ways  of  obtaining  good  stabilization  are  studied. 

One  by  use  of  a  gas  rocket;  which  yields  a  non  linear  system  ("On  -Off"  type) 
studied  by  the  phase  plane  method.  The  other  by  use  of  flywheels:  the  stabiliza¬ 
tion  is  then  a  linear  one.  In  both  cases  it  is  shown  that  a  great  simplification  can 
be  introduced  in  the  stability  equations  by  the  introduction  of  some  approximations. 
Even,  in  the  non  linear  case,  the  approximations  (practically  always  justified)  brin 
an  enormous  simplification  in  the  computations  by  supressing  the  coupling  between 
the  three  axis  around  which  the  stabilization  is  applied. 
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CHAPTER  I 
Introduction 


The  goal  of  this  thesis  is  the  study  of  the  attitude  stabilization  of  an  earth 
observation  satellite*  The  orbit  is  assumed  circular  or  near  circular  (eccentricity 
e  <  10’*);  and  the  altitude  is  above  200  miles,  in  which  case  the  atmospheric  drag  is 
negligable*  Also*  the  perturbation  torques  are  assumed  very  small  (and  indeed  they 
are)*  compared  to  the  inertial  torques  on  the  satellite  -  and  in  doing  so  it  will  be 
possible  to  linearize  the  equations  of  stability  (small  angles  approximations:  cos  8 1 , 
sin  6  ~  6).  The  case  of  the  insertion  into  orbit  where  large  errors  appear  will  not  be 
studied. 

After  a  first  part,  devoted  to  the  study  of  the  orbit  from  a  mechanical  point  of 
view,  the  main  part  of  this  report  will  discuss  the  stabilization  of  the  satellite  by  use 
of  rockets  (Chapter  III). 

The  corrective  torques  needed  for  stabilization  and  produced  by  the  rockets  will 
be  of  the  *  on-off*  type,  with  deadzone  and  time  delay  (in  practice  this  is  a  good  approx¬ 
imation  for  most  cases). 

The  study  of  such  a  system  in  three  dimensional  space  is  rather  complicated. 

But  by  the  use  of  suitable  approximations,  it  's  possible  to  make  the  study  in  one 
dimensional  space. 

Then  some  particularities  of  the  stabilization  by  a  non-linear  system  will  be 
studied  "reticence  behaviour*  phenomena,  limit  cycle.  It  will  appear  that  generally 
two  orientations  modes  (coarse  orientation  mode,  and  fine  orientation  mode)  are  needed 
and  can  be  obtained  by  switching  the  gain  of  an  amplifier.  Also,  a  quick  study  of  the 
amount  of  stored  gas  for  the  rockets  is  made,  for  a  satellite  life  of  8  to  10  months. 

As  such  stabilization  by  rockets  is  not  precise  enough  if  the  satellite  carries 
cameras,  a  following  chapter  will  be  devoted  to  a  flywheel  stabilization  system. 

To  conclude  it  can  be  said  that  both  systems  are  needed:  the  rocket  stabilization 
for  large  perturbations;  and  the  fly  wheel  stabilization  to  achieve  the  stabilization  with¬ 
in  the  required  limits. 
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CHAPTER  II 
Mechanics  of  the  Orbit 

2.1  Syetem  of  Reference* 


System  of  Reference 
Fig.  2.  1 


a)  X'f  Y' ,  Z'  are  fixed-body  axes  in  the  satellite.  For  simplification  sake  they 
are  chosen  to  coincide  with  the  principal  axces  of  inertia  of  the  vehicle.  The  reduction 
in  generality  of  the  results,  in  so  doing,  is  not  serious,  since  most  of  the  satellites  have 
some  sort  of  symmetryjand  control  will  be  desired  about  these  principal  axes. 

The  vehicle  is  built  and  oriented  in  such  a  way  that  (X*,  Y* ,  Z' )  almost  coincide 
with  the  reference  system  XYZ  which  we  shall  describe  now. 

b)  If  X,  Y,  Z  (see  Fig.  2. 1)  are  chosen  to  be  almost  parallel  to  X' ,  Y',  Z'  then  the 
variation  angles  will  be  small,  and  we  shall  be  able  to  linearize  equations  by  the  following 
approximations 

sin  0  ~  9 
cos  9  ~  1 


*  See  Reference  (1) 
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So,  as  shown  in  Fig.  2.1,  the  Z  axis  is  aligned  to  the  geocentric  vertical  and  taken 
positive  outward  from  the  center  of  the  earth.  The  X  axis  is  in  the  plane  containing 
Z  and  the  velocity  vector,  and  oreinted  in  the  direction  of  the  velocity  (for  a  perfectly 
circular  orbit  X  and  the  velocity  vector  will  coincide)  Y  is  chosen  to  form  a  right-hand 
system  with  X.  Z. 

Such  a  reference  system  is  not  fixed  in  space,  but  rotates  with  an  angular  vel¬ 
ocity  f). 

For  a  circular  orbit  we  get 

[2  =  £1  =  0 
x  z 

O  =  D  =  constant 

y  O 

2.  2  Disturbances  Torques 

Disturbance  torques  can  be  put  into  two  categories,  internal  and  external. 

2.2.1  Internal  disturbances  are  due  to  the  motion  of  some  internal  piece,  a  shutter  a 
rotating  magnetic  recorder, etc.  and  will  not  be  taken  into  account  for  generally  they 
are  very  small. 

2.2.2  External  disturbances  can  be  divided  into  secular  and  cyclic  types: 

Secular  disturbances 

a)  Aerodynamic  drag  which  is  of  the  order  of  1  dyne  cm  at  the  altitude  considered 
(at  200  miles).  Above  200  miles  aerodynamic  forces  are  negligable,  at  least  for  a 
satellite  life  of  8  to  10  months. 

_  3 

b)  Meteoritic  collision  for  a  me* ec rite  of  more  than  10  gm,  the  probability  of 

# 

occurence  is  quite  small  . 

c)  Earth's  magnetic  field 

The  torques  are  due  to  the  eddy  current  effects  or  permanent  magnet  effects  in 
the  satellite  which  react  with  the  earth  field.  These  may  be  one  of  the  most  important 
sources  of  disturbances.  They  may  be  minimized  by  some  special  precautions.  For- 
instance  in  the  case  of  eddy  currents,  by  using  a  skin  of  high  resistivity  or  a  skin  made 
of  small  sections  not  electrically  connected,  also,  internal  current  loops  can  be  made  by 
pairs  of  twisted  wires  whenever  possible.  The  effects  of  the  permanent  magnets  will 
be  corrected  by  compensating  coils  along  the  3  orthogonal  axes.  The  order  of  magnitude 
of  the  magnetic  torques  is  100  dyne  cm  to  200  dyne  cm. 

d)  Differential  gravity 

If  a  principal  moment  of  enertial  is  not  aligned  with  the  yaw  axis,  a  torque  T^  is 

$ 

See  Reference  (2). 
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created.  For  a  circular  orbit  and  a  small  variation 


where 


Tgx  =  AI  3w‘ 


AI  =  I -  I 

e)  Earth's  oblateness.  ’ 

The  oblateness  of  the  earth  is  a  perturbation  which  acts  directly  on  the  satellite 
orbit.  The  plane  of  the  orbit  oscillates  slowly,  the  amplitude  depending  upon  the  in¬ 
clination  ij/of  the  orbit  upon  the  equator.  But  such  a  perturbation  is  generally  negligible. 


Cyclic  perturbations 

a)  Solar  radiations 

5  2, 

The  solar  radiation  pressure  is  of  the  order  of  5  x  10  dynes/  cm  . 

b)  Differential  gravity. 

Differential  gravity  effects  the  solar  cells  paddles. 

c)  Magnetic  fields 

For  all  of  the  above  cyclic  disturbances  the  frequency  is  about  twice  the  orbital 
frequency  (the  period  of  the  orbit  being  around  100  minutes  for  an  altitude  of  200  miles) 
and  the  maximum  torque  is  some  scores  of  dyne,  cm  (100  to  200  dyne,  cm  is  a  good 
average  value).  Finally  only  the  action  of  the  gravity  gradient  is  simple  enough  to  be 
accounted  for  simply  in  the  equations  of  motion.  The  other  perturbations  being  in¬ 
dependant  of  the  oreintation  of  the  satellite,  will  be  considered  as  those  unavoidable 
disturbances  which  occurs  in  every  serovomechanism.  Their  action  will  be  that  of  an 
impulse  (meteoritic  impact),  a  step  (magnetic  disturbances),  or  a  sinusoid  (cyclic 
perturbations). 


In  any  case  they  are  shown  to  be  small  as  far  as  their  effects  are  concerned; 
and  the  perturbations  on  the  stabilization  will  be  well  within  the  tolerance.  For  instance, 
a  perturbation  of  200  dyne  x  cm  (which  is  high)  will  give  the  satellite  an  angular 
acceleration  of 


e  = 


(I  >108  g.cm.2) 


or  V  «.  10‘®  rad. /sec. 2 

which  is  negligible  for  the  period  of  time  we  shall  consider  (we  shall  see  that  the  oscil¬ 
lations  have  a  period  of  about  100  sec.  ). 

lft 

2.3.  Equations  of  motion 

Two  different  sets  of  equations  can  be  developed  depending  on  how  the  correction 


♦See  Reference  (1)  and  (2) 
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torques  are  produced;  by  ejection  of  a  certain  mass  of  gas,  or  by  an  internal  motion 
of  a  flywheel.  If  there  are  no  moving  parts  in  the  satellite  then  the  motion  of  the  vehicle 
in  orbit  will  be  described  by: 


and 


H  =  [i]  w 

Zi  =  e  +  q 


(2.1) 


Where  0:  angular  velocity  of  the  reference  system 

6:  angular  velocity  of  the  body  w.  r.t.  the  reference  system 
u:  angular  velocity  of  the  body  w.r.t.  the  inertial  space 
w.  r.t.  :  an  abreviation  for  "  with  respect  tb"  . 


0 

0 

I 

z 


Now  if  we  call  T  any  torque  (disturbance  or  control  acting  on  the  vehicle)  we  have  the 
relation 


dH 
cTF~  - 


T 


The  expansion  of  equation  (2.1)  and  (2.2)  is: 


ie  +  40  (I  -  I  )0  +  O  (I  -  I  +  I  )  0  a  T  +  T  , 
xx  o'y  z'x  o'x  y  z'  z  xc  xd 


I  6  +  30  (I  -  I  )  0 
y  y  o  x  z  y 


-  Tyc+Tyd 


10  +0  (1-1)8  -O  (I  -1+1)9  =  T  +  T  , 

z  z  o  y  x'  z  o'  x  y  z'  x  zc  zd 


(2.2) 


(2.3) 


A  detailed  derivation  of  these  equations  is  given  in  Reference  1.  These  equations  are 

valid  for  a  circular  (or  near  circular)  orbit  where  O^  =  0q  =  constant 

0=0=0 
x  o 

Txc  is  the  control  torque  along  x  axis.  Tx(j  is  the  disturbance  torque  along  x  axis, 
(note  that  the  gravity -gradient  torque  is  included  in  the  left-hand  side  of  the  equations). 
If  there  are  moving  parts  in  the  satellite  (flywheels)  then: 

H  .  =  H  ,+  E(H  .  ,  ) 

syst.  veh.  wheels' 
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or  .more  explicitly: 


H.y.r[K)'"v.h.  4-'"v.h. 


and 


H  =  T 
syst. 


where 


H  syst 
[K] 


=  Angular  momentum  of  the  whole  system. 

=  Inertia  of  the  satellite  minus  flywheels. 
=  Inertia  of  a  flywheel. 


(2.4) 


The  expansion  of  (2.4)  and  (2.5)  is  given  in  references  (see  Ref.  (3)).  For  the  resulting 
equations  see  Chapter  IV. 
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CHAPTER  III 


Non-Linear  Stabilization  Using  Jets 

We  shall  study  the  stability  problem  of  a  satellite  using  jets  to  create  a  con¬ 
trol  torque.  The  jets  operate  on-off,  and  the  control  system  is  nonlinear. 

3.1  Nonlinear  Stabilization  System  in  the  Three  Dimensional  Case. 

a)  Stability-  In  the  absence  of  control  torque  the  equations  (2.3)  after  using  La¬ 
place  transformed  yield: 


[lxS  +4«0  V  h\ 

l#x+no«V  V*I»)fl-" 

A(s) 

(3.1) 

f  I  S2+  3 ft  2<I  -  I  ) 

]e 

B(s) 

(3.2) 

l  y  O  X  z ' . 

J  y 

\  I ,S2+fl  2  (I  -  I  ) 
l  z  o  y  x' 

1  8  -  Q  S(I  -I  +  I  )8  = 

J  z  o  x  y  z'  x 

C(s) 

(3.3) 

where  A(s),  B(s),  C(s)  stand  for  the  Laplace  transform  of  the  perturbation  torques  and 

the  initital  conditions  (8  (o),  8  (o) . ).  The  system  will  be  stable  if  the  roots  of  the 

x  y 

following  equations  don't  have  a  positive  real  part: 

Vs'*  “.X-  V  ■ 0 

is2* «„ V i.' y y  . 0 

-syvyy  i/.c/iy  y 

or 

IyS2  +30,2  (Ix-  Iz)  =  0 

i  i  s4+s2(i  «  2(i  -  i  )+4i,n  2(i  -  i„) 

X  Z  X  o  y  X  c  O  y  z 

+s2n  2  (i  -  i  + 1  )2  +4fl  4  (i  -  i  )(i  -  i  ) 
o'x  y  z '  o  y  z  y  x 

An  evident  solution  is  that 


I  -  I  >0 
x  z 


I  -  I  >0 
y  x 


I  -  I  >0 
y  x 


or 

I  >1  >1  (3.4) 

y  x  z 

Such  a  relation  is  sufficient  but  not  necessary. 


b)  Approximations ;  and  reduction  to  the  one  dimension  case 

Assuming  now  that  the  inequality  (3.4)  is  satisfied  assuring  that  the  system  is 
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stables  the  details  of  stability  can  be  investigated.  For  example,  from  equation  (3.  3) 
we  notice  that  the  damping  or  rate  term  is  non-existant;  and  a  critical  study  of  equations 
(3.  1)  to  (3.  3)  will  show  that  no  damping  effects  are  present.  For  example,  the  mag¬ 
nitude  of  the  terms  in  Eq.  3.1  will  be  examined.  The  value  of  1  ,  for  a  compact  vehicle 

8  9  2  * 

of  several  hundred  pounds,  is  of  the  order  of  10  ,  10  g  cm  (for  a  compact  satellite, 

which  is  generally  the  case,  the  variation  of  I  from  one  principal  axis  to  the  other  is  no 

8  8*' 

more  than  1  order  of  magnitude  -So,  I  ~  10  ,  1^  ~  10  ,  is  determined  from:  L=d.T= 

I  0  where  T  is  the  thrust  of  the  control  jets  and  d  the  distance  between  the  jets.  Usually 

x  x  -  3  14'* 

T  is  of  the  order  of  some~10  dynes.  V/ith  d=100cm  one  obtains  j—  ~  10  ,  and  0  = 

10  rd/  S^.  Also  -  10  ^  assuming  a  rotation  of  the  satellite  around  the  earth 

in  100  minutes.  The  quantity  (I  -  I  )  is  of  the  order  of  I  ,  at  the  most  as  2  or  3x10  rd  and 
2  A  y  ^  ^  o  x  * 

4U  (1  -  I  )  0  ~  10  1  at  the  most.  Also  §<  10  '  rd/  s,  and  in  average  0  —10"  rd/  s 

o  y  z  x  x  '  6  ' 

(such  values  can  be  obtained  from  any  articles  in  the  litterature  giving  data  of  the  per¬ 
formances  cf  any  satellite  ).  Lastly,  1^(1^-  1  +  *z)®x  <  10"^  1^. 

•  • 

It  appears  that  all  terms  of  equation  (3.1)  are  negligible  compared  to  Ix®x*  The 
same  for  (3.2)  and  (3.  3).  These  equations  can  now  be  approximately  written  as: 

I  9  =  A(s)  +  T  (s) 

x  x  x' 

10  =  B(s)+Ty(s)  (3.5) 

10  =  C(s)  +T  (s) 

Z  * 

T(s)  is  the  Laplace  transform  of  the  control  torques.  The  coupling  can  be  neglected  and 
each  equation  can  be  studied  separately.  Now,  there  is  no  rule  which  will  allow  one  to 
choose  T(s),  but  as  it  will  be  seen  later,  for  the  linear  system,  the  form  of  T(s)  needed 
will  be:  Tx(s)  =  A9x+  B  0  .  Thus,  a  good  first  try  for  the  nonlinear  system  is:  aTx(s)  & 
sign  (A©x+  B©x)  or  better  Tx(s)  =  sign  (A©x+  ®®x+  C),C  being  a  deadzone  which  always 
appears  in  such  systems.  Nonlinear  control  and  compensation  is  discussed  in  detail  in 
the  next  section. 

3.  2  Study  of  a  Nonlinear  System  in  the  One  Dimension  Case  (by  the  phase  plane  method) 

The  block  diagram  of  such  a  nonlinear  gas  jet  system  can  be  represented  in  the 
following  way  (See  Fig.  3. 1).  perturbation 


Fig.  3.1:  Block  diagram  of  a  jet  stabilization  system  (stabilization  around  one  axis  only) 
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The  parameter  2a  is  the  deadzone  width  which  exists  in  most  on-off  systems.  T  is 
the  sum  of  different  delays  which  may  occur  in  the  loop  between  the  detection  of  a 
deviation  and  the  applied  control  correction.  Kp  is  the  gain  of  the  rate  g^jro.  K  is  the 
gain  of  the  sun  and  earth  sensor  which  detects  any  deviation  6.  We  call  -jj? *  X.  Now 
as  we  have  a  stabilization  system,  the  input  will  be  zero  (a  correction  of  attitude  com¬ 
manded  from  the  earth  can  be  considered  as  a  perturbation),  and  the  loop  can  be  con¬ 
sidered  in  the  following  way  (Fig.  32). 

PERTURBATION 


_-TS 


uU 


IS 


xs+ 1 


L*-| — 


-Q 


Fig.  3.2  Loop  diagram  for  jet  stabilization  around 
one  axis 


To  begin  with,  the  perturbations  will  not  be  considered  (they  will  be  introduced  later), 
a)  t  =  0 

The  open  loop  can  be  translated  into  the  equation: 


‘  *  L  sign  (-0  -  X  tt-  +  a) 

dr  dt  “ 


(3.6) 


The  phase-plane  method  will  be  used  to  study  the  system.  Two  new  variables  are  in¬ 
troduced  in  (3.  6) 


Then  (3.6)  becomes: 


x  =  -9 

,  dx  d0 

y  =  +3F  =  '3F 


I  ^  =  L  sign  (x+  Xy  +  o) 


(3.7) 


Sucha  differential  equation  can  be  represented  by  a  curve  (C)  in  a  plane  of  coordinates 
(x,  y)  called  the  phase  plane.  The  representative  point  (x,y)  of  Eq.  37  will  move  on  C. 

When  the  sign  of  x+Xy  +  a  will  change  there  will  be  a  discontinuity  of  shape  on 
C  (whose  differential  equation  will  leap,  for  instance,  from^  1  to  ^  *  -1). 

Such  discontinuities  occur  on  the  lines  x+Xy  +  a*  0,  which  are  called  the  switching  lines. 
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The  equation  of  curve  (C)  in  the  phase  plane  will  now  be  examined.  Assuming,  at  the 
start,  that  sign  (x+\y-a)<0,  Eq.  (3.  8)  becomes  I  =  -L. 


or  since 


dx  dy  dy  dx  dy 

y=-at’3F  =  cfx  '  HF  =  3x  y 


I  y 


dy  _ 

as 


=  -  L 


which  yields  by  integration 


Iy  =  -  2 Lx  +  Constant 


Assuming  that  at  t=  0,  x=xq,  y=0then:  Iy  =  -2L(x-xq) 


(3.8) 

(parabola ). 

(3.9) 


and  the  time  to  coast  along  such  a  parabola  is  t=  -  -j-  y  which  is  obtained  by  replacing 
dx=ydt  in  (3.8).  Assuming  equation  (3.9) 


L  Jl  ,,  . 

T  t  =  '2(x-  xo} 


b)  T  ^  0 


Now,  due  to  the  delay  r,  there  is  no  switch  on  the  line  y=  ^  .  But  the 

point  (x,  y)  will  move  on  the  parabola  (3.9)  during  t  more  seconds,  before  the  appli¬ 
cation  of  the  torque  takes  effect  (see  Fig.  3.4). 


Fig.  3.4  Influence  of  a  time  delay  t 
on  the  switching  line 
(Due  to  the  delay  t  the  switching  occurs  at  (2)  instead  of  (1)). 


From  the  following  equations : 


t,=  - 


r  yl 

21  ,  . 

■it  (xr xo) 


xi+  Xy^-  Q=  o 


We  find 


11 


1 

tj*  -X  +(X2+^-  (xq-  a))*,  -X  +A 


(the  other  value  of  tj  is  <  0  and  so  cannot  be  considered). 
Now  we  know  that  tj+T  or  t^  =  -X  +  A  +  t  . 

Substituting  the  above  value  of  t2  in  the  following  equations 


h~ 


I 

L  y2 


*  21  .  . 
*2=  '  T7 (x2*  xo> 


yields  the  relation  between  y^,  and  given  in  Eq.  3.  10. 

Lt 


x 2-  a  -  7T  (t  -2\  ) 
(t-  X) 


(3.10) 


which  is  the  equation  of  a  delayed  switching  line.  Eq.  (3.10)  is  valid  for  a  point  (x.  y) 
being  in  the  fourth  quadrant  (or  the  first).  Symmetrically,  i.e.  ,  in  changing  a  into  -a 
•y-into — y  (Fig.  3.5  results) 


We  find  the  switching  line  for  the  second  quadrant 

Lt 

xt  a  t  TT  (t  -  2X ) 


(T-  X) 


(3.11) 
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After  switching  on  the  lines  (3*10)  or  (3.11)  the  jets  will  be  off.  i.e.  no  control  torque 

d0 

will  be  exerted  and  y  (or  ^-)  will  be  constant.  In  the  resulting  motion  the  describing 
point  (x,  y)  will  move  on  a  line  parallel  to  the  x  axis,  (see  Fig.  3.6). 

Due  to  the  delay  r ,  the  switch  jets  on  command  will  occur  r  seconds  after  point 
(3)  in  the  phase  plane  is  reached. 


We  have 

and 

from  which  we  deduce  y4= 


X4  a  Xj+  T  y4 


y4s  v 


Xj+  Q 


x4  -  a 


which  is  the  equation  of  the  third  switching  line. 


Changing  a  into-a  in  this  exp/es^lon  gives  the  4th  switching  line. 
Finally  we  have  4  delayed  switching  lines 


(3.12) 

(3.13) 


i-.-O 
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The  Fig.  (3.6)  makes  clear  that  the  slope  of  the  switching  lines  must  be  <  0,  if  we 
don't  want  the  describing  point  to  diverge,  i. e. ,  an  unstatie  system  to  occur.  There 
fore,  stability  requires  X  >0  since  r  is  always  «  X. 

3.  3  Particularites  of  a  Jet  Stabilization  System  -  (The  single  dimentional  case) 

A  study  of  Fig.  3.  6  shows  two  particularities  of  such  a  system, 
a)  "Reticence*  behavior  (the  use  of  such  a  term  will  become  clear  later). 


Fig.  3.7  Reticence  behavior 


Fig.  3.8  Reticence  behavior  -  its  causes 

As  the  trajectory  coasts  through  the  deadzone  (see  Fig.  3.7)  and^reaches  a 
switching  line,  a  negative  torque  is  commanded  (position  (1)  Fig.  3.7);  *  -y  de¬ 

creases  but  8  is  still  increasing  due  to  the  inertia  of  the  satellite  -  and  if  X  is  small 
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the  sign  of  0  +  X  0  +  a  may  change  before  0  changes  in  sign.  The  rockets  will  then  be 

—  JQ 

switched  off  with  a  unchanged  in  sign  (position  (2)).  So  the  same  pair  of  jets  will 
again  be  switched  on  (position  (3))  after  a  certain  increase  of  0  etc.  In  short,  the  same 
pair  of  jets  changes  the  direction  of  rotation  of  the  satellite  after  many  on  cycles  and 
not  after  one  unique  big  thrust. 

This  way  of  stabilizing  the  satellite  may  seem  economical,  for  instead  of 
oscillating  around  the  equilibrium  position  the  satellite  comes  very  near  to  it,  very 
quickly.  It  is  important  to  study,  however,  what  the  conditions  must  be  imposed  on 
the  control  system,  in  such  a  case. 


Fig.  3.9  Limit  cycle  behavior 


The  graph  (3.8)  shows  that  at  some  point  the  |  slope  |  of  the  parabola  must  become  great¬ 
er  than  the  j  slope  |  -i-  of  the  switching  line.  This  point  must  be  in  the  range  of  values  of 
-2 

y<  2x10  rd/s  (a  greater  value  can  come  only  from  a  bad  insertion  -  accident  which  is 

_3 

not  taken  in  account)  and  y  <  10  rd/s  (below  such  rotation  no  control  is  needed:  for 
such  small  rotation  is  not  a  trouble  for  the  observation  instruments  of  the  satellite.) 

Finally:  j  slope)  of  parabole  =  j-  -i 


and  we  want 


or 


1 

X 


< 


L  1 

Ty 


The  lowest  value  of  y  is  around  10 


X  ~ 


10 


15 


•o 


L 

T 


>  10 


-4 


But  in  most  cases,  ^  ~  1  or  2  xlO”^  rd  (sec  . 

So  the  reticence  behavior  will  occur  only  for  very  small  values  of  y  when,  as  we 
shall  see  later,  the  jet  must  give  way  to  flywheels. 

b)  Limit  cycle  behavior 

Fig.  3.7  shows  that  the  point  (x,  y)  is  becoming  "nearer  and  nearer"  to  the 
origin  but  due  to  the  symmetry  of  the  switching  lines  it  may  fall  into  what  is  called  a 
"limit  cycle"  ,  i.  e. ,  the  point  (x,  y)  will  endlessly  follow  the  same  trajectory  around  the 
origin.  See  Fig.  (3. 9). 

The  coordinates  of  points  (1)  and  (II)  are: 


or 


xl=  Xlf  Ua+a+7T  <T‘  2X>1 

Xra 

Y1~'YIIS  :  “  ymax 

T  -  X 

Y  s  +  Lt  (t~  ZX) 

max  "  51  (t-X) 


(3.14) 


(3.15) 


X  ^  is  such  that 
max 

y2  =  -  (X  -  X  ) 
’  I  '  max' 

in  which  x  and  y  are  replaced  by  x^  and  y^  to  find 
The  result  is 


Ym«  ‘it 


+  X  =  a+£r 
—  max  4T 


+  P(t-X) 


with  p  «r  T(T.,--2X) 
t-  X. 


Note  that  t 
X 


10  sec 

5  to  15  sec  »  t 


(3.16) 


so  P  *2t 
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Fig.  3.10  Two  ways  of  approaching  the  limit  cycle  (1)  staying  "outside"  of  it  and 
(2)  Moving  inside  and  outside  the  limit  cycle. 

A  graphical  study  (see  Fig.  3.10)  shows  that  the  trajectary  approaches  the  limit 
cycle  in  the  steady  state.  It  may  sometimes  "enter*  the  zone  delineated  by  the  limit 
cycle  curve,  but  it  will  get  out  quickly  and  start  again  approaching  such  limit  cycle. 

3.4  Study  of  a  Perturbation  on  the  Rocket  Stabilization  System 

This  paragraph  will  be  devoted  to  the  study  of  the  damping  of  x  and  y  versus 
time,  i.e.  ,  the  finding  of  a  relation  between  x  (and  y)  at  the  points  (0),  (6)  on  Fig. 

3.  11  and  the  time.  (The  variation  of  x  (or  y)  between  points  (0)  and  (6)  can  be 
considered  as  the  variation  driving  a  psuedo  period). 

Indeed  the  perturbation  started  at  a  point  (x.,  y.)  somewhere  in  the  plane;  but 
it's  more  practical  to  start  from  a  particular  point  for  instance  from  a  point  on  a  switch¬ 
ing  line);  and  to  go  back  to  (x^  y{)  from  the  chosen  particular  point  by  a  change  in  the 
time  origin.  The  following  computations  will  give  the  values  of  the  coordinates  of  points 
(o),  (1),  (2);  from  there  it  will  be  possible  to  deduce  easily  3,  4,  5  by  use  of  the  sym- 


17 


Fig.  3.11  Phase  Plane  for  x(* -9)  and  y(» ) 


metry  in  the  Fig.  3. 11.  The  following  equation*  will  be  used: 

parabola  y^»  +  (X-X0)  (3.19) 

switching  line  y*  +  )  (3.20) 

which  can  be  approximated  by  y*  — since  r  is  generally  one  hundred  times  less  than  X , 
and  since  the  4  switching  lines  can  be  approximated  by  2  switching  lines  only.  The  time 
along  a  parabolic  path  varies  as  given  by  the  equation 

t«  |y|  (3.21) 


The  intersection  of  a  switching  line  (which  replaces  two  switching  lines  almost 
in  coincidence)  and  of  a  parbola  yields  the  points  (0)  and  (2),  or  from  Eq.  (3.19)  and 
(3. 20)  it  turns  out: 


A 


V 


LX 


+  A 


(X  is  known  to  be  negative) 
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with  A*  [(“j~  )2 

^Note  that  A  is  smaller  than 
Also  (3.  22)  yields 


1 

(a-x^f 

for  Xj  greater  than  a 


y 


_  LX. 

r  t 


+  A 


] 


Form  Eq.  (3.21)  the  time  to  move  from  point  (0)  to  point  (2)  is  t--  t  such  that  t,-t  =* 


The  final  result  is: 

y0+y2  = 


(3.22) 


+X  (ly0l+ly2l> 

2  LA. 


(3.23) 

(3.23bis) 


From  there  the  graph  of  y  can  be  drawn.  For  if  y  is  known  (and  it  is  known  from 
the  initial  perturbation  y^)  Eq.  (3.23)  gives  y2  and  then  (3.23)  gives  the  time  elapsed  be¬ 
tween  the  values  y  and  y, ;  and  (3.  21)  shows  that  y  is  linear  in  t.  From  there  point  (3) 

O  b 

is  deduced: 


since 


y3=  y2 


t,  -  t„ 


— 

y2 


dx 

dt 


constant  =  y^ 


a  integration  yields:  ^  variation  of  x  _  2a. 

y2  y2 

But  the  quantity  t3-  t2  is  negligable  compared  to  t2*  tQ  (a  is  of  the  same  order  of 
magnitude  as  y2);  so  the  2  points  (2)  and  (3)  will  be  in  coincidence  when  drawing  the 
graph  oi  y  versus  time.  As  for  the  points  (3),  (4),  (5)  they  are  deduced  from  the  pre¬ 
ceding  computations  merely  by  replacing  a  by  -a,  y  by  -y-  and  the  point  (0),  (1),  (2) 
respectively  by  (5),  (4),  (3)  i.e.  ,  it  turns  out  that: 

y,+  v^  -2-^  (3. 24) 


VS’T  <ly5l+|y3l) 


and  y*  —  t, 


from  Eq. 


(3.  21),  i.e. ,  y  ir  varying  linearity  with  respect  to  time,  and  ine 
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•lope  of  y  versus  t  is  y  . 

From  point  (5)  to  (6)  y  remained  constant  but  the  time  t,  -  t-  is  negligable  (being 

2a  ”  3 

equal  to  — ). 


Finally  it  appears  that  the  graph  of  y  versus  time  will  be  made  of  segments  of 
lines,  of  slopes  +  y  ,  and  the  extremities  of  such  segments  will  be  given  by  Eq.  (3.22) 
or  (3. 24).  See  Fig.  3. 1 3.  Note  that  the  units  in  Fig.  3. 13  has  been  normalized,  i.e., 
instead  of  y  and  t,  y=  y.^-^  and  ^  are  used;  the  slopes  of  die  lines  becom  then  +  1,  and 
the  difference  between  2  maxima  of  y  is  4  1;  and  the  time  elapsed  between  Z  maxima  of 
y  is  equal  to  the  sum  of  the  absolute  values  of  these  maxima.  It  must  be  noted  that  when  a 
positive  maximum  of  Y  is  less  than  1,  then  the  following  maximum  of  Y  should  be  pos¬ 
itive  in  order  to  satisfy  Eq.  3.22;  that  is  in  contradicition  with  what  is  shown  in  Fig.  3.11 
where  2  consecutive  maximum  of  y  (and  so  of  Y)  are  of  opposite  signs.  This  contradiction 
is  due  to  the  fact  that  the  approximation,  made  in  assuming  2  switching  lines  in  coincidence 
is  no  longer  valid  as  it  appears  in  the  following  Fig.  3.  12. 

Case  one  the  2  switching  lines  (a)  and  (b)  can  be  assumed  in  coincidence.  Case  two  such  an 
approximation  is  no  longer  valid. 


The  integration  of  y*  yields  x.  In  fact  from:  Y*  gjr*  2X^T  *  ^ ‘  T~ 


the  integration  gives 


4LX 

~T 


2  -T 

_/  YdT 
o 


and  since  Y  is  linear  in  T,  f  Y  dT  is  the  equation  of  a  parabola.  Being  careful  about 
the  concavity  of  the  parabola,  it's  easy  to  draw  the  graph  of  x  versus  T.  (See  Fig.  3.14). 


Fig*  3.14  Damping  of  the  Satellite 

Angular  Error  Vereue  Time 
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The  two  figures  of  x,  and  y  versus  t  have  been  drawn  one  under  the  other,  because  it  is 
easy  to  sketch  one  graph,  knowing  the  other.  Note  that,  exactly,  the  portions  of  parabola 
of  Fig  3.14  should  be  recorded  by  straight  lines  originated  by  the  Integration  of  y*  constant 
(See  Fig.  3.13);  but  is  has  been  shown  that  such  constant  values  of  Y  are  maintained  dur¬ 
ing  too  short  a  time  to  be  taken  into  account. 

It  must  be  noted  also  that  th/?  parabolas  are  symetrical  around  the  Y  axis;  and 
that  they  are  connected  between  themselves  on  the  T  axis  (this  is  shown  by  Fig.  2.11, 
where  the  change  of  sign  occurs  on  the  horizontal  lines,  between  2  parabolas.  The 
maxima  of  X  are  computed  from  the  graph  of  Y  versus  T. 

The  2  curves  3.13  and  3.14  give  some  conditions  on  the  stability  parameters.  For, 
to  have  a  quick  damping  of  a  perturbation: 

X  should  be  large  (so  for  a  given  T,  t  will  be  small) 
should  be  as  small  as  possible. 

The  preceding  study  was  made  to  get  an  idea  of  the  damping  of  6  and  6  when  the 
perturbation  is  an  impulse.  The  system  is  stable,  and  provided  a  good  choice  of  the  par¬ 
ameter  -j-,  X  etc. . .  the  damping  will  be  fast  enough.  For  other  types  of  perturbations 
(step,  sinusoid)  the  phase  plane  method  is  not  very  adequate.  Nevertheless,  if  the  per¬ 
turbations  are  not  too  big  (i.e. ,  some  100  dyne,  cm)  they  can  be  compared  to  an  impulse: 
for  under  such  disturbances  (a  step  6  and  0  will  be  built  up  very  slowly;  and  beyond  a 
certain  deadzone  the  stabilization  system  will  act  swiftly  enough  so  that  it' s  possible  to 
neglect  the  small  influence  of  the  step  disturbance  on  6  and  6  . 


3.  5  Energetical  Study 

When  the  angular  velocity  is  varying  from  0  to  y,  the  energy  used  is 

1  ,  2 

2 


Forinstance,  in  the  case  of  the  Fig.  3.13  it  will  be: 

' 2  '  ‘/5 


E  =7Iy02  +Iy22  +Iy52+Iy72 


or  since  Y  *  y  • 


ZlA 


-i-2  E=[(^Yo2+(Y0-1)2  +  (Y0-2)2+Y0-3)2  ] 


e  =  ri 


<IYoZ-llYo+17> 
7  2  llyo<2LX> 

ryo - : — 


17  <i^)2] 


The  energy  used  for  a  given  perturbation  will  be  smaller  if:  -£•  is  large  and  X  small. 
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CHAPTER  IV 
Flywheel  Stabilization 


4. 1  Equationa  of  Motion 

Referring  to  Chapter  II  and  Reference  3  the  following  equationa  are  derived  for 
the  motion  of  the  vehicle  ami  its  flywheels 


ixs2+  (y  2+  jry  \  J nzs |dx-  y y «  s-  j flys 
-  J(fizs+  nx» )  |  i  s2  j"  j{nxs-  nzu ) 

-i  -  i +i  )  w  s+ jn  s  -jn  si  i  s2 -hi  - 1  )  w2+ jn  « 

xyz'  y  I  xiz  y  x'  y 


W 


+  J 


o 

s 

0 


u 

0 

s 


["]■  [*] 


(4.1) 


With:  Ix>  I  ,  Iz  inertia  of  the  satellite  minus  flywheels  on  the  X,  Y  and  Z  axis. 

J  inertia  of  each  flywheel. 

ft  angular  velocity  of  the  satellite  in  the  X,  Y,  Z  frame. 

to  angular  velocity  of  the  X.  Y,  Z  frame  with  respect  to  an  inertial  frame. 
T  disturbaces  corques. 


Since  the  gravitational  torques  are  known,  T  can  be  partitioned  into  two  parts: 


with  T 


grav. 


T  _  t  +  T 

dist.  grav. 

(I  -  I-)w2  0  0 

y  2  2 

0  (I  -  I  Xo  0 
'  X  z' 

e  oo 


In  a  very  rough  approximation  of  equations  (4-1)  it  is  possible  to  neglect  almost  all  of 
the  left-hand  side  terms  of  these  equations.  For  instance,  the  first  equation  of  (4-1)  can 
be  written  as : 


e  + 

XX  x 


[tty-  Iz)  4w  2  +  jny«]  +  0y  JO, 

+  ®2[(Ix-  Xy+  V10  -  J°y]  +J(V“  nz> 


(4.2) 


=  T 


dist. 
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*  Reasonably"  it  can  be  assumed: 

.  JO  «  I*  (as  J«  10‘5  I  it  would  imply  fl  «  103  rd|S  which  is  justified) 

.  6  <  10  rd|S  all  greater  values  being  reduced  by  the  rocket  stabilisation. 

•  •  .5  ,2 

.  6  >10  rdJS  which  is  true  only  in  certain  cases. 

.  u  «  10'3  rd  |s 

Then  Eq.  (4.2)  can  be  rewritten  as: 


i  V  +  jh 

XX.  X 


diet. 


and  the  same  for  the  2  other  coordinates,  J ft  is  the  control  torque;  it  can  be  chosen  as 


J  f2  a  K  6  +  K  6 
P 

such  choice  being  justified  by  the  fact  that  it  will  appear  convenient. 


4.2  Study  of  a  Flywheel  Stabilization  around  One  Axis. 

With  such  approximations  the  coupling  between  axes  disappears.  From  now  on 
the  stabilization  around  each  axis  can  be  studied  separately.  The  block  diagram  of  a 
stabilization  by  flywheel  will  be  the  same  for  each  axis.  It  has  been  drawn  in  Fig.  4.1; 


T  =  DISTURBANCE 


Fig.  4.1  Block  diagram  of  a  stabilization  by 
flywheel  system 


«  meaning  *  far  less  than". 
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0  :  correction  on  the  attitude  satellite  commanded  from  earth  (in  general  6  «  0). 
c  c 

is  the  correction  on  the  orientation  of  the  Z  axis  (see  Fig.  2.1).  It's  due  to  the 

rotation  of  the  satellite  frame  around  the  earth;  and  in  case  of  a  circular  orbit  0  .  * 

TCI  e 

ut.  t  being  the  time,  and  "w  "  a  constant.  The  transfer  function  of  the  flywheel  has  been 
taken:  |^r-  •  The  flywheel  being  mainly  an  electric  motor,  this  is  the  transfer  function 
of  some  electric  motor.  But  the  justification  is  given  in  Appendix  4,  and  I  is  the  in¬ 
ertia  of  the  satellite  (minus  flywheel)  around  the  axis  (I  stands  for  I  or  I  or  I  ).  The 

x  y  z 

transfer  function  of  diagram  (4-1)  is  easily  found  to  be 


6  a 


1+  TS 


(1+TS)  +  A(KpS+K) 


Td+^ 


A(K  S  +  K) 


IS  (1+TS)  +  A(KpS+K) 


ref. 


a)  Tj=  0  (no  disturbance),  as  the  orbit  is  almost  a  circle: 

de 


ref 

"3T~ 


=  w  *  Constant 


we  want  0  -  ®re£  a  0  for  a  ramp  input  of  0  ^ 


£  =  0  ref. -9 


v  _  lim  c  A 

jCVy"”  g  ^  o  ^  •  — -y 


-  .  (K+K  S)  oo 
JSfc(l+TS)  p 


i,  e. ,  £  —  0 

So  no  error  is  due  to  the  steady  variation  of  6  with  respect  to  an  inertial  frame, 
b)  Influence  of  disturbances  T^ 

ea  -  it !§. -  .  T.  (4.4) 
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The  stability  of  the  system  is  studied  by  the  Routh  criterion.  The  characteristic 
equation  of  0  is 
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The  system  will  be  stable  if  -TK  >  0,  i.  e.  >  T.  (In  practical  case  =  \  is  of 
10  to  20).  The  study  of  (4-4)  is  difficult  enough  due  to  the  3rd  degree  denominator. 

If  yet  it's  assumed  that  the  delay ypjfj  is  small  enough  i.e.,  T  =  0  we  get  a  greatly 
simplified  expression: 


8 


T 


d 


The  study  is  then  straight-forward,  and  can  be  found  in  any  feedback  litterature  (see  for- 
instance  Ref.  3).  Such  a  study  can  only  give  an  idea  of  the  values  for  KjKp,  T  etc. :  for 
the  approximations  to  get  (4.4)  for  instance  are  not  always  justified.  In  particular,  in 
(4.4)  there  is  no  mention  of  the  coupling  of  the  flywheels;  and  yet  it's  evident  that  such 
wheels  behave  as  gyroscopes.  H.  Cannon  in  the  ARS  Journal  (See  Ref.  3)  has  studied  an 
interesting  way  of  coupling  the  three  flywheels  so  that  their  momenta  can  be  exchanged  as 
the  satellite  rotated  around  the  earth  and  so  a  great  amount  of  energy  can  be  spared. 

As  the  Routh  criterium  has  shown  that  the  system  was  stable  (if  >  T .  K) ,  the 
*  flywheel  stabilization  will  be  useful  for  small  and  high  frequency  perturbations  (where 
the  gas  rocket  may  be  destabilized). 

Both  systems  will  be  needed  for  a  good  stabilization;  the  flywheel  system  furnishings 
the  damping  of  the  perturbations  already  damped  by  the  gas  rocket  system. 


26 


CHAPTER  V 
Conclusion 

In  conclusion  it  can  be  said  that  in  order  to  achieve  a  good  stabilization  of  the 
satellite,  two  corrective  systems  are  needed.  One.  in  most  cases  a  gas  rocket  stabil¬ 
ization  system,  to  correct  strong  perturbations  torques  ("strong"  ,  meaning  here:  of 
the  order  of  some  hundreds  dyne  x  cm). 

Such  a  system  is  also  very  interesting  for  impulse  disturbances;  and  the  phase 
plane  method  is  very  adequate  to  study  the  system  under  such  perturbations.  In  fact, 
in  the  study  done  here,  all  the  disturbances  have  been  reduced  to  impulse  disturbances  by 
assuming  that  step  disturbances  and  sinusoidally  varying  disturbances  are  small  and  of 
large  periods;  then  they  build  up  9  and  9  very  slowly  (due  also  to  the  large  inertial  of 
the  satellite);  and  the  gas  rocket  system  will  act  only  when  9  and  0  (when  the  satellite  was 
only  under  the  influence  of  disturbances  torques.)  Such  a  method  -  reducing  all  distur¬ 
bance  torques  to  an  impulse  is  only  valid  for  slowly  varying  and  small  perturbations.  It's 
justified  for  most  of  the  perturbations  encountered  in  space.  But  an  interesting  study 
could  be  done  on  the  influence  of  a  high  frequency  perturbation  on  the  gas  rocket  system, 
or,  in  other  words,  under  what  conditions  may  the  satellite  be  destabilized;  then,  the 
method  of  the  first  harmonic  should  probably  be  used  instead  of  the  phase  plane  method. 
Also  for  low-level  disturbances  (for  instance  some  dyne  x  cm)  the  coupling  between  the 
3  axes  cannot  be  neglected;  and  the  study  of  the  gas  rocket  system,  with  coupling  be¬ 
tween  the  3  axes  of  stabilization,  is  quite  impossible  without  any  mechanical  means  of 
computation.  Also  low-level  disturbances  need  not  be  studied  due  to  the  fact  that  a  gas 
rocket  system  is  ineficient  for  the  correction  of  a  low  level  disturbance  because  of  its 
limit  cycle,  therefore,  a  flywheel  stabilization  system  has  been  studied  for  small  distur¬ 
bances. 

In  that  case  the  system  is  a  linear  one;  but  the  coupling  between  the  3  axes  of 
stabilization  is  difficult  to  neglect  and  must  be  thoroughly  examined.  The  flywheel 
system  should  be  studied  by  means  of  a  computer.  Otherwise  the  equations,  though  not 
impossible  to  solve,  become  rather  difficult.  An  interesting  feature  of  a  flywheel  system 
is  that  a  coupling  between  the  3  flywheels  allow  a  great  economy  of  power;  for  the  energy 
can  be  transferred  from  one  flywheel  to  another  without  too  much  losses.  Such  coupling 
has  been  studied  by  H.  Cannon  in  the  ARS  Journal  (see  Ref.  3).  For,  if  the  complete 
study  of  the  stabilization  of  a  satellite  is  difficult,  many  separated  studies  can  be  entirely 
worked  and  without  too  much  trouble,  even  without  the  use  of  a  computer. 
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APPENDIX 


4.  Derivation  of  the  transfer  function  for  a  flywheel 
A  diagram  of  such  wheel  is: 
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Fig.  4.  I  A  Flywheel  System 


Assuming  0  =  Q 
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it  turns  out  the  following  relation: 
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Equations  (2)  and  (1)  then  yield: 
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